December 14, 2009 



CASIMIR EFFECTS NEAR THE BIG RIP SINGULARITY IN VISCOUS 

COSMOLOGY 

Iver Brevik 1 *, Olesya Gorbunova 2 ^, Diego Saez-Gomez 3 * 
1 Department of Energy and Process Engineering, 
Norwegian University of Science and Technology, N-7491 Trondheim, Norway 
2 Tomsk State Pedagogical University, Tomsk, Russia and 
3 Institut de Ciencies de I'Espai (IEEC-CSIC), E-08193 Bellaterra (Barcelona), Spain 

Analytical properties of the scalar expansion in the cosmic fluid are investigated, especially near 
the future singularity, when the fluid possesses a constant bulk viscosity f. In addition, we assume 
that there is a Casimir-induced term in the fluid's energy-momentum tensor, in such a way that the 
Casimir contributions to the energy density and pressure are both proportional to I /a 4 , a being the 
scale factor. A series expansion is worked out for the scalar expansion under the condition that the 
Casimir influence is small. Close to the Big Rip singularity the Casimir term has however to fade 
away and we obtain the same singular behavior for the scalar expansion, the scale factor, and the 
O ' energy density, as in the Casimir-free viscous case. 
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I. INTRODUCTION 

From a hydrodynamical viewpoint it is almost surprising to notice that the cosmic fluid - whether considered in 
the early or in the late epochs - is usually taken to be nonviscous. After all, there are two viscosity (shear and bulk) 
coefficients naturally occurring in general linear hydrodynamics; the linear approximation meaning physically that 
one is considering only first order deviations from thermal equilibrium. The shear viscosity coefficient is evidently 
of importance when dealing with flow near solid surfaces, but it can be crucial also under boundary-free conditions 
such as in isotropic turbulence. In later years it has become more common to take into account viscosity properties 
of the cosmic fluid, however. Because of assumed spatial isotropy in the fluid the shear viscosity is usually left out; 
any anisotropic deviations like those encountered in the Kasner universe become rather quickly smoothed out. Thus 
qq , only the bulk viscosity coefficient, called £, remains in the energy-momentum tensor of the fluid. One should here 
note, however, that at least in the plasma region in the early universe the value of the shear viscosity as derived from 
kinetic theory is greater than the bulk viscosity by many orders of magnitude. Cf., for instance, Refs. [3 and 0. 

Early treatises on viscous cosmology are given by Padmanabhan and Chitre [|[ and Gr0n [3| , the latter paper being 
an extensive review of the field up to 1990. We have ourselves dealt earlier with viscous entropy production in the 
early universe Q, and viscous fluids on the Randall-Sundrum branes [U, Hj]. Cataldo et al. considered viscous dark 
energy and phantom evolution using Eckhart's theory of irreversible thermodynamics pj- Kofinas et al. considered 
the crossing of the w = — 1 barrier using a brane-bulk energy exchange model with an induced gravity curvature 
correction Q. As discussed by Nojiri and Odintsov Q and by Capozziello et al. [nj, a dark fluid with a time 
dependent bulk viscosity can be considered as a fluid with an inhomogeneous equation of state. Some other recent 
papers on viscous cosmology are Refs. (ill Il2l [l3j . 

A special branch of viscous cosmology is to investigate how the bulk viscosity can influence the future singularity, 
commonly called the Big Rip, when the fluid is in the phantom state correspondin g to the thermodynamic parameter 
w being less than —1. Some recent papers in this direction are Refs. [hj, |l5|, [l6|, [l7|, [l8|, llEj, 12CJ, |2l|. In particular, 
as first pointed out in Ref. [14], the presence of a bulk viscosity proportional to the scalar expansion 9 can cause the 
fluid to pass from the quintessence region into the phantom region and thereby inevitably lead to a future singularity. 

The purpose of the present paper is to generalize these viscous cosmology theories in the sense that we take into 
account the Casimir effect. We shall model the Casimir influence by writing the total Casimir energy in the same 
form as for a perfectly conducting shell, identifying the cosmological "shell" radius essentially with the cosmic scale 
factor. This is a very simple, though natural, approach to the problem. The approach is of the same kind as that 
followed in an earlier quantum cosmology paper dealing with the expanding FRW universe in the nonviscous case 
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[22j . We mention that there are also other ways of treating the influence from the Casimir effect in cosmology; the 
reader may consult, for instance, Refs. [23l [24l |25|. [2rl [27| . 

As bulk viscosity corresponds to inhomogeneous Hubble rate dependent terms in the effective equation of state, the 
situation is quite similar to that of modified gravity theory. For a general treatise on modified gravity, see the review 
of Nojiri and Odintsov [281 ]. 



II. FORMALISM, WHEN THE CASIMIR EFFECT IS OMITTED 

We start with the standard FRW metric, 

ds 2 = -dt 2 + a 2 (i)dx 2 , 



(1) 



and set the spatial curvature k, as well as the cosmological constant A, equal to zero. We let subscript 'in' refer 
to present time quantities, and choose t in = 0. The scale factor a(t) is normalized such that a(0) = a in = 1. The 
equation of state is taken as 



P = wp, 



(2) 



with w constant. As mentioned, w < — 1 in the phantom region (ordinary matter is not included in the model). The 
bulk viscosity £ is taken to be constant. The energy- momentum tensor of the fluid can be expressed as 



= pUfjJJv + p(g l/tv + U^Uy) 



(3) 



where £/ M is the comoving four- velocity and p = p — (9 is the effective pressure, 9 = 'iH = 3d /a being the scalar 
expansion. The Friedmann equations take the form 



The energy conservation equation leads to 



9 2 = 24irGp, 

9 + 1.9 2 = -4nGp. 
6 

P +(p + p)9 = C9 2 . 



(4) 
(5) 

(6) 



From the above equations the differential equation for the scalar expansion, with w and £ as free parameters, can be 
derived as 



9 + -(1 + w)9 2 - 12t:GC9 = 0. 
The solution is, when subscript zero signifies that the Casimir effect is so far left out, 



ot/tc 



l + i(l+™)0 in i c (e*/tc_i) ' 



Here 9i n is the initial (present-time) scalar expansion and t c the 'viscosity time' 

1 

tr 



(7) 



(8) 



(9) 



12ttGC 

Since (f + w) < by assumption, it follows from Eq. ((8]) that the future singularity occurs at a rip time t s o given by 

2 1 



tso — tc In 



1 - 



This means that in the nonviscous limit C —> 0, 



ts 



1 + w 9 in t c 



1 



1 + W 9 m ' 

which is independent of £. At the other extreme, in the high viscosity limit t c — > 0, 

-2 f 
1 + w 9 in t c 

showing that t S Q becomes small. The fluid is quickly driven into the Big Rip singularity. 



t s o — t c In 



(10) 



(ff) 



(12) 
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III. THE CASIMIR EFFECT INCLUDED 

As mentioned above, a simple and natural way of dealing with the Casimir effect in cosmology is to relate it to the 
single length parameter in the (k = 0) theory, namely the scale factor a. It means effectively that we should put the 
Casimir energy E c inversely proportional to a. This is in accordance with the basic property of the Casimir energy, 
viz. that it is a measure of the stress in the region interior to the "shell" as compared with the unstressed region on 
the outside. The effect is evidently largest in the beginning of the universe's evolution, when a is small. At late times, 
when a — > oo, the Casimir influence should be expected to fade away. As we have chosen a nondimensional, we shall 
introduce an auxiliary length L in the formalism. Thus we adopt in model in which 

where C is a nondimensional constant. This is the same form as encountered for the case of a perfectly conducting 
shell 29]. In the last-mentioned case, C was found to have the value 

C = 0.09235. (14) 

The expression ([13]) is of the same form as adopted in Ref. [22;] (cf. also [3(3]). It is strongly related to the assumptions 
made by Verlinde when dealing with the holographic principle in the universe [3l|. Cf. also the papers [32j and [33| 
dealing with the holographic principle applied to viscous cosmology. 

In the following we shall for definiteness assume C to be positive, corresponding to a repulsive Casimir force, though 
C will not necessarily be required to have the value (TT4"]) . The repulsiveness is a characteristic feature of conducting 
shell Casimir theory, following from electromagnetic field theory under the assumption that dispersive short-range 
effects are left out (|29|; cf. also [I(|). Another assumption that we shall make, is that C is small compared with 
unity. This is physically reasonable, in view of the conventional feebleness of the Casimir force. 

The expression (|13p corresponds to a Casimir pressure 

, - - 1 dE < C (15) 

Pc ATr(La) 2 d{La) AttLW V ' 

and leads consequently to a Casimir energy density p c oc 1/a . 
The Casimir energy-momentum tensor 

T; v = PcU^Uv + p c ( 9fiV + U^U„), (16) 
together with the Casimir equation of state p c = w c p c , yield the energy balance 

— + (1 + w c )9 = 0, (17) 

Pc 

having the solution p c a z ^ 1JrWa ^ = constant. To get p c oc 1/a 4 we must have w c = 1/3. The Casimir contributions to 
the pressure and energy density become accordingly 

C 3C 
Pc= 8^' Pc= MV- (18) 

The Friedmann equations now become 

i + ^ = -4^(„„-C. + i ^), (20) 

while the energy conservation equation preserves its form, 

p+(l + w)p8 = (6 2 . (21) 



Note again that we are considering the dark energy fluid only, with density p and thermodynamical parameter w. 
The ordinary matter fluid (dust) is left out. 
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Solving p from Eq. (|19p and inserting into Eq. (|21[) we obtain as governing equation for the scalar expansion 

1 SGf 1 
6+-(l + w)9 2 -l2irG(9 = -(l-3w)^ [ . (22) 

It is convenient to introduce the constant a, defined as 

a = -(l + w)>0, (23) 

and also to define the quantity X(t). 

X(t) = 1 - iaf? m i c (e t/tc - 1), (24) 

which satisfies 

X(0) = 1, X(t s0 ) = 0. (25) 
In view of the assumed smallness of C we now make a Stokes expansion for 9 to the first order, 

0(t) = 9 {t) + C9 1 (t)+O(C 2 ), (26) 
using henceforth the convention that subscript zero refers to the C = case. The zeroth order solution is 

9 (t)=e in X- 1 e t ^, (27) 
in accordance with Eq. ([8]). It corresponds to the zeroth order scale factor 

ao(t)=X-&, (28) 

satisfying ao(0) = 1 as before. 

As the right hand side of Eq. (|22]) is already of order C, we can replace a(t) with ao(t) in the denominator. Thus 
we get the following equation for the first order correction coefficient 0\\ 

3G 



(aB in X- 1 e t l u + 12^GC) 9i = -(1 - 3w)^jX&. (29) 



We impose the same initial condition for the scalar expansion as in the C = case: 9(0) = 0q(0) = 9i n . It means 
according to Eq. fl2B) that 6>i(0) = 0. 

The homogeneous solution of Eq. (|29p . called 6\h, may be written 



0ih(t) = exp 

satisfying ^(O) = 1. The the full solution becomes 



/ (a9 m X- 1 e t/t ^ + 12irG()dt 
Jo 



(30) 



3G „ , , f'll 



(t) = -(l-3w)— [ 9i h (t)- I —dt, (31) 



°lh 



satisfying ^i(O) = 0. The two terms on the right hand side of Eq. (|26|) . 9o(t) and C9\(t), are accordingly determined. 



A. Analytic approximation for low viscosity 

Although in general the expression for 9(t) has to be calculated numerically, the main features of the solution can be 
shown already analytically. Consistent with the assumed smallness of C we need not distinguish between the rip time 
t s corresponding to C ^ and the rip time t s o corresponding to C = 0. Let us assume for mathematical simplicity 
the low- viscosity limit 



9% n t c S> 1, 



(32) 
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being physically the most important case also. It corresponds to t s o/t c — 2/(a6i n t c ) -C 1. Then, 

X(t)~la6 in (t s0 -t) = t -?^, (33) 

O (t)~--^— , (34) 
a r s0 - t 

»o«)~(^)*. (35. 
Both 6*o (i) and do(i) diverge (a > by assumption). Using Eq. (|3"3"|) we can calculate 9\(t) from Eq. ([3T]) . 

W V ; 2L 4 8 + 9a (1-^Aao) 2 

From the expansion (|26p we thus obtain for the scalar expansion to the first order in C, 

' W - ^T^W 1 ( 3 } 4^ 8T9^ } • ( 3? ) 

The viscosity is absent in this expression. This is as we would expect in view of the low-viscosity approximation. 

The expression (|3T|) cannot, however, be valid near the singularity. The reason is that the Taylor expansion in C 
in Eq. (|26p is not applicable at t = t S Q . The solution (f3"T|) can be applied safely as long as t stays considerably smaller 
than t s o. By making a first order expansion in t/t s o of the expression between the curly parentheses we can write the 
solution in simplified form as 

% 2 1 f ,..„•> 3GCat s n 1 ,„„. 
^^^T^l 1 -^- 3 -)^^*}. *<*- (38) 

As (1 — 3u>) > this means that 9(t) becomes slightly reduced by the Casimir term. The repulsive Casimir force 
causes the energy density p in Eq. (fT9")l to be slightly smaller than in the C = case. 

To deal with the conditions close to the singularity, we have to go back to the governing equations themselves. 

B. Behavior close to the singularity 

To begin with, it is instructive to list the general classification of possible future singularities as given in Refs. [34| 
and [2l[. If t s denotes the rip time, one has four types, 

(i) Type I (original "Big Rip"): For t —>t s , a — ► oo, p — > oo, and \p\ — » oo, or p and p are finite at t = t s . 

(ii) Type II ("sudden"): For t — > t s , a — > a s , p — > p s , and |p| — » oo, 

(iii) Type III: For t — ► t S) a — > a s , p — > oo, and |p| — > oo, 

(iv) Type IV: For t — > t s , a — > <z s , p — > 0, |p| — > 0, or p and p are finite. Higher order derivatives of H diverge. 

The singularities we have been contemplating above are seen to be of Type I. As a — > oo near the singularity, we 
can draw the important conclusion from Friedmann's equations (1191) and (|20[) that the influence from the Casimir 
term fades away. Close to the singularity we simply obtain the same solutions as in the Casimir-free case. 

Let now t^ s denote the singularity time in the presence of viscosity. We thus get 

t sC =t c ln(l + --L). (39) 
It corresponds to 9(t s () = oo. We see that t S £ is always less than the singularity time for the nonviscous case, 

t si < i Slf =o = (40) 



For the scalar expansion we find close to the singularity, again assuming for simplicity low viscosity so that 9i n t c 1 
0. 

»(*) = 7^T, (41) 



In turn, this corresponds to 



C. On the nonviscous case 

It may finally be worthwhile to consider the entirely nonviscous case, while keeping C > 0. Setting £ = we get 
from Eq. (|22[) the governing equation for the scale factor a: 

1 GC 
a 3 a+ -(1 + 3w)a 2 a 2 = -(1 - 3w) — . (44) 
^ zL 



This equation is still not solvable analytically. We get for X(t), 6o(t) and ao(t) the same expressions as in Eqs. (j33|) - 
([35]) . Similarly, we get for 8{t) the same expansion ([37| as before, with t s0 = 2/(a0i n ). These results are as expected, 
in view of the property of continuity with respect to variation of parameters. 

IV. CONCLUDING REMARKS 

Considerable attention has recently been devoted to the behavior of the dark energy fluid near the future singularity. 
Various possibilities have been contemplated. In addition to the references given above [23| -[27|, we may refer also 
to the papers 13511 . It has even been suggested that the future singularity can be avoided via quantum gravity effects. 
Thus in Ref. [361 ] it is shown how the universe may turn into a de Sitter phase. 

Let us finally summarize our results above: 

1) If £ > and C = 0, the rip singularity time t s o is given by Eq. (|10p. In particular, if £ — > then Eq. (fTT|) holds. 

2) If £ > and C > 0, the scalar expansion 9(t) is given by the first-order series (|26p . with 9o(t) and 6%(t) given 
by Eqs. (f2"T)) and (|3"Tj) . In the low-viscosity case 9i n t c ^> 1, 9(t) is given by the series (|3"T|) which, however, is not 
applicable near the singularity as 9(t) is not analytic in C at the singularity. 

3) Near the singularity, the Casimir effect fades away and the viscous rip time t s ^ is given by Eq. ([55)) . Corresponding 
values for 9(t),a(i) and p(t) near the singularity follow from Eqs. (|4"Tj) - (T4"51) . 

4) If £ = and C > 0, the governing equation for a(t) is Eq. (JUJ). 
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